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1 .  INTRODUCTION 


•tlJ 


Recent  advances  in  sem i conductor  growth  techniques  have  led  to  a  wide¬ 
spread  interest  in  the  physics  of  ultra-small  semiconducting  systems.  Quantum 
wells,  superlatt ices ,  double-barr ier  resonant  tunneling  structures,  and  a 
variety  of  other  exotic  structures  have  become  the  objects  of  extensive  inves¬ 
tigation  [1,2].  The  interest  in  these  ultra-small  systems  is  motivated  by  two 
factors.  First,  their  optical  and  electrical  properties  have  quasi-two- 
dimensional  (2-d)  features  which  frequently  offer  distinct  advantages  in 
device  applications.  Second,  actual  physical  systems  whoso  electron  dynamics 
are  quasi-2-d  provide  one  with  a  rich  testing  ground  for  theoretical  models. 


One  of  the  systems  attracting  considerable  attention  is  the  doubl e-barr ier 
resonant  tunneling  (DBRT)  structure  [3].  A  typical  structure  consists  of  two 
thin  -50  A  AlGaAs  layers,  separated  by  an  equally  thin  (-50  A)  GaAs  layer,  all 
of  which  are  embedded  in  a  single  GaAs  crystal .  The  regions  to  the  left  and 
right  of  the  barriers  (usually  beyond  spacer  layers)  are  n-doped  and  usually 
electrically  contacted  for  transport  studies.  Current-voltage  characteristics 
of  this  device  show  an  enhancement  in  the  current  when  the  appl  ied  voltage 
aligns  the  quasi-Fermi  energy  of  incoming  electrons  with  the  energies  of  the 
quasi-bound  states  in  the  quantum  well  region.  A  number  of  theoretical  calcu¬ 
lations  have  been  done  to  describe  the  nonlinear  current  response  [*J-l6]  in 
this  system.  However,  controversy  still  exists  regarding  the  basic  mechanism 
behind  the  nonlinearity  [1*1,15,17-20].  Despite  the  large  number  of  studies  to 
date,  we  have  not  found  in  the  literature  any  calculation  of  the  local  density 
of  states  (DOS)  for  even  a  highly  simplified  model  of  a  DBRT  structure.  The 
local  DOS  provides  information  about  resonant  states  and  gives  one  a  quantita¬ 
tive  measure  of  the  extent  to  which  the  dynamics  are  quasi-2-d.  The  purpose  of 
this  paper  is  to  present  such  a  calculation. 


In  section  2,  using  a  simple  model  potential,  we  calculate  the  eigenvalues 
and  eigenfunctions  of  an  effective  mass  Schroedinger  equation.  These  results 
are  used  in  section  3  to  compute  the  local  DOS  in  the  quantum  well  region  of 
the  potential.  In  section  A,  we  relate  this  local  DOS  to  an  integral  of  a 
one-dimensional  ( 1  — d )  DOS.  This  latter  DOS  has  sharp  peaks  at  energies  corre¬ 
sponding  to  the  quasi-bound  states  between  barriers. 


2.  THE  MODEL 


We  consider  a  simplified  model,  which  is  defined  by  the  following  effec¬ 
tive  mass  Hamiltonian: 


V*  +  V(z)  , 


where  mn  is  the  effective  mass  of  electrons  at  the  bottom  of  the  GaAs  conduc¬ 
tion  band  and  the  double  barrier  potential,  V(z),  is 


V(z)  =  VQ.5(z  *  a )  *  VQ<5(z  -  a) 


In  this  model,  the  two  AlGaAs  potential  barriers  have  been  replaced  by  delta 
function  barriers  of  strength  V  ,  separated  by  a  distance  2a  along  the  z-axis 


(the  growth  direction).  The  parameter  VQ  is  given  by 


VQ  =  bAVc  , 


where  AV  is  the  conduction  band  discontinuity  and  b  is  the  barrier  width. 


We  solve  the  one-electron  Schroedinger  equation 


HY(r)  =  Ef(r ) 


subject  to  periodic  boundary  conditions  in  the  x-  and  y-directions 


f(x+L,  y ,  z)  =  vCx,  y,  z) 


Y(x,  y+L,  z)  =  f(x,  y ,  z) 


Since  equation  (4)  is  separable,  we  can  write  the  wavefunction  in  the  product 
form 


f(r)  =  -  e 


—  P 


ip(  z )  , 


where  k^  =  (kx,ky,o),  kx  =  27rnx/L,  ky  =  2irny/L,  and  nx,  ny  take  the  integer 
values  0,  ±1,  ±2,....  The  z-part  of  the  wavefunction,  \|j(z),  satisfies  the 


reduced  equation 


r  '(z)  *  iTJ—  Ce  -  v(z)]^(z)  =  0  , 


where  e  =  E  -  fi2k?/2m  .  For  ij>(z),  we  choose  the  vanishing  boundary  conditions 


♦(-  ?)  =  o  . 


!);(-)  =  0  . 


The  presence  of  the  delta  functions  imposes  two  jump  discontinuity  condi¬ 
tions  on  the  derivative  of  i|i(z)  at  the  delta  function  positions,  which  are 
given  by 


ij/'(-a  +  0 +  )  -  \|)'(-a  -  0  +  )  =  Y<ji(-a)  , 


i|;'(a  +  0  +  )  -  ip’  (a  -  0  +  )  =  Y«p(a )  , 


_ “SS 


■V.  yOvUVik  ■A.kC’A  k m  kv  k.  kv  k«  k.  k.  kv  k.k  .sk.  a  k» kv_k  .ok.  k.k 


where  Y  =  2mcV0/-h2,  and  0+  is  a  positive  inf initesimal  quantity.  These  two 
conditions  may  be  found  by  integrating  equation  (7)  over  the  infinitesimal 
intervals  (-a  -  0+,  -a  +  0+)  and  (a  -  0  +  ,  a  +  0+).  In  addition,  we  require 
the  wavefunction  to  be  continuous  at  the  delta  function  positions 


( -a  -  0  +  )  =  ip  ( -a  +  0  +  )  , 

4>(a  -  0  +  )  =  »p( a  +  0  +  ) 


(10a) 


(10b) 


To  look  for  a  solution  to  equation  (7),  we  take  advantage  of  symmetry  and 
solve  the  eigenvalue  problem  in  the  region  0  <  z  <  L/2.  We  then  look  for  even 
wavef unct ions  of  the  form 


\  A. (k)cos(kz) 


0  <  z  <  a 


*  u(z)  =  i  L 

/  A~(k)cos(kz)  +  A,(k)sin(kz)  ,  a  <  z  <  - 


and  odd  wavef unctions  of  the  form 


^  B  (k)sin (kz) 


0  <  z  <  a 


<|»„(  2)  -■  ,  U2) 
°  (  (k  )cos (kz )  +  B^(k)sin(kz)  ,  a  <  z  <  -  . 

For  the  even  wavef unctions ,  when  we  impose  the  conditions  in  equations  (8b), 
(9b),  and  (10b),  we  get  a  set  of  three  homogeneous  equations  for  A 1 ( k ) ,  A2(k), 
and  An(k).  To  have  a  nontrivial  solution,  we  require  the  determinant  of  the 
coefficient  matrix  to  vanish.  This  gives  an  equation  for  the  allowed  k's 
which  label  the  even-parity  states: 


cos(ka)sin(k£,  -  ka)  +  cos(kil)  =  0 


(13a) 


Here  we  have  used  the  convenient  definition  H  =  L/2.  Solving  the  three  homoge¬ 
neous  equations  for  the  ratios,  we  find 


A3(k)  Y  2f,  ^ 

I^kT  "  k  cos  (ka)  ’ 


(13b) 


A  (k) 

— y — r-  =  1  -  —  sin(ka)cos(ka) 
A^ (k )  k 


(13c) 


Applying  the  same  boundary  conditions  to  the  odd  solutions  given  in  equa¬ 
tion  (12)  and  setting  the  determinant  of  the  coefficient  matrix  to  zero,  we 
obtain  the  equation  satisfied  by  the  allowed  wavevectors  labelling  the  odd- 
parity  states: 


-  sin(ka)sin(kl  -  ka)  +  sin(k!,)  =  0 
k 


(19a) 


Solving  the  associated  system  of  equations  leads  to 


B2(K)  t  2 

5700  ■  ■  k  sln  <ka>  ■ 


(14b) 


B3(k)  Y 

g ~^k ^  =  1  +  sin(ka)cos(ka) 

The  constants  A,  and  B,  are  evaluated  from  the  normalization  condition 


(14c) 


/: 


*ctk 


(z)  | 2  dz  =  1  , 


where  a  =  e  (o)  for  the  even  (odd)  solutions.  We  find 


1 


A?  ( k ) 


=  l{ 


1  + 


sin(2ka ) 
2k  2, 


+  (]  ~  !)[&) 2sin2(2ka)  -  i  sin  ( 2ka )  ♦  Q'cosMka)] 


1 

kl 


cos  k(2+a)sin  kti-a)^-^-^  sin2(2ka)  -  ^  sin(2ka)  -  cos“(ka)  +  ij 


~  ^  [ik  sin(2ka)  “  ijcos2 (ka )sin  k(2+a)sin  k(2-a)j 


(15) 


and 


B2  ( k ) 


=  11 


sin(2ka) 

2k2 


(l  ~  l)  i<  |sin(2ka)  +  ~  [sin4(ka)  +  ^  sin2(2ka)]| 


+  cos  k(2+a)sin  k(2-a)  [^sin‘*(ka)  -  ^  sin2(2ka)J  -  ~  sin(2ka)  -  1 

~  Vl  [|]<  sin(2ka^  +  ^  sin2  (ka  )s  in  k(2+a)sin  k(2-a)|  .  (16) 


Th°  energy  oigonvajues  associated  with  the  eigenfunctions  in  equation  (6)  arQ 
given  by 


ka 


2m 


(17) 


With  vanishing  boundary  conditions  on  il(z)  at  z  =  ±2,  the  limit  of  zero 
strength  delta  functions  (Y  -*■  0)  is  identical  to  the  limit  where  the  delta 
functions  are  piaced  on  the  boundaries  (a  ■*  2).  In  both  cases  one  recovers 
the  1-d  particle  in  a  box  problem,  where  A,/A1  =  Bp/B1  =  0,  A0/A1  =  B^/B,  =  1, 
and  A,  -  B,  =  (2/L)l/2.  J 


3.  THREE-DIMENSIONAL  DENSITY  OF  STATES 


We  now  construct  the  single-particle  Green's  function 


c(?.?'!E)  -  111  -2S - 2* — r 


:*  a  k  E-E^+iO 
k|  a  k 


where  the  eigenstates  y  +  (r)  are  given  by  equations  (6),  (11),  and  (12;, 

CiK 

=  (kj^,  ka)*  Here  a  (equal  to  e  or  o)  labels  a  state's  parity,  ki 

k  ),  and  ka  are  given  by  the  roots  of  equations  (13a)  and  ( 1  ha )  for  even- 
odd-parity  states,  respectively.  Using  the  Green's  function,  we  ea.cuiat. 
local  DOS  (including  both  spins): 


D(z,E)  =  -  -  Im  G  ( r  ,  r ,-  E ) 


U  l  l  l  Uok  (z)|2  «(E  -  E-) 
r*  a  k  a 

kj_  a 


In  the  limit  when  the  system  size  goes  to  infinity,  with  "a"  he;-, 
stant,  the  density  of  allowed  wavevectors  becomes  2tt/L.  This  a.iows  t 
change  the  sums  in  equation  (20)  to  integrals: 


D  ( z  ,  E )  = 


ll  TT 


i:  -i  i 


dk  l  |*ah(z)|>  «(E  -  E-;  . 

0  a 


The  integration  over  k^  =  (k^.k^)  is  over  all  positive  and  negative  wav.. 

tors,  whereas  the  z-component  wavevector,  k,  is  integrated  over  pos : 
values  only.  Using  the  explicit  form  of  the  wavefunctions  and  changin 
spherical  momentum  coordinates,  we  have  for  the  region  -a  <  z  <  +a 


D ( z , E) 


_ 1 _ 

®aoEo  Jo 


/2)( E/E0) 1 


/  2 


dq  [Fe(q)cos2^  q)  +  FQ(q)sin2(^  q) 


where 


Fe  ^ C1  ^  q2  +  U2cos2(q)  -  Uq  sin(2q) 


I'  *55 


Fo(q)  =  q2  *  U2 


sinT(q)  +  Uq  sin(2q) 


i  iU  f, 

'  *  25  -  e 


i2q- 


f  18) 


a  n  d 


•  ve  - 
t  .  ve 

’  y  .  - 

.?  ?a  ) 

23b) 


9 


S? 


the  dimensionless  potential  strength  U  is  defined  by  U  =  Ya  = 
convenient  energy  scale,  EQ  =  ir2ti2/8mca2 ,  has  been  introduced 
F  (q)  and  F  (q)  are  related  to  the  wavefunction  amplitudes  by 


2mVQa/fi  ,  and  a 
The  functions 


lim  iA2fa\  =  F  (q) 

1/1*0  Ua'  e 

lim  =  F  (q) 

.  1  W  o  M 


a/2,+0 


The  local  DOS  given  by  equation  (22)  is  a  sum  of  two  terms.  The  term  contain¬ 
ing  Fe(q)  gives  the  local  density  of  even  states  in  the  well,  and  the  tern 
with  FQ(q)  gives  the  local  density  of  odd  states.  Each  function  F  (q),  ct  =  e 
or  o,  is  an  even  function  with  an  infinite  sequence  of  (complex  conjugate) 
pairs  of  simple  poles  in  the  complex  q-plane  (see  fig.  1).  Each  pair  of  poles 
with  Re ( q )  >  0  corresponds  to  one  resonance.  In  the  limit  U  -*■  0  (no  delta- 
function  barriers),  these  poles  move  away  from  the  real-q  axis  to  infinity  ana 
Fa(q^  -  1  on  the  real  axis.  In  this  limit,  D(z,E)  is  simply  the  local  DOS  for 
a  free  electron  gas  in  a  box  of  volume  2aL2 . 


—  0.5 
-1.0 

-1.5 

-2.0 


5.0  10.0  15.0 

Re  q 


20.0 


Figure  1.  Locations  of  pol'es  of  Fg(q)  are  shown  for  U  ^ 

3.  Poles  of  F  (q)  (not  shown)  lie  between  poles  of  Fg(q) 


As  we  approach  the  limit  of  strong  barriers,  U  >>  1  ,  the  poles  of  the 
functions  Fa(q)  move  in  toward  the  real  axis  in  pairs  (one  pole  above  and  one 
pole  below)  (see  fig.  1).  In  the  region  of  U  >>  1,  and  0  <  q  <  U,  the  func¬ 
tions  Fa(q)  are  well  represented  by  a  sum  of  Lorentzians 

F  (q)  -  l  - - ^ -  ,  (24) 

n  (q  -  q  ) 2  +  T2 
an'  an 

where  the  real  and  imaginary  positions  of  the  poles,  qan  ±  iTan,  are  functions 
of  U.  For  U  -*■  «  we  have  ran  ->•  0,  qgn  >  (2n  +  1  ) tt/2 ,  and  qQ“  -*•  mr.  In  this 
limit  we  find 


lim  F  (q) 

ii  a 

U*°° 


-  I  6(<I  -  qan} 

n 


Rather  than  look  at  the  local  DOS  in  more  detail,  we  consider  its  integral 
over  the  well  volume 


r/2)(E/E0] 1/: 


N(E)  = 


KT2fl2  J  0 


dq  Fe(q)  ♦  F0(q)  +  Fe(q)  -  F0(q) 


in(2q) 

2q 


The  function  N(E)  gives  the  number  of  states  in  the  well  per  unit  energy 
interval.  In  the  limit  of  weak  barriers,  U  <<  1 ,  we  find 


U-0  1  / 2m\ 3/2 


2aL2/E  , 


which  is  the  DOS  f or  a  free  electron  gas  in  a  volume  2aL2. 

In  the  limit  of  strong  barriers,  U  >>  1 ,  we  find 

N(E)  ~ ^  i  e ( E  -  n2Eo)  ,  (27) 

n  =  1 

where  9(x)  =  1  for  x  >  0  and  0  ( x )  =  0  for  x  <  0.  This  is  the  well-known 

staircase-like  DOS  one  would  expect  in  a  quasi-2-d  system.  For  intermediate 
values  of  U,  the  DOS  N(E)  is  plotted  in  figure  2. 


Figure  2.  Dimensionless  density  of  states  between  barriers, 
(u^fiVmL2 )n( E)  ,  is  plotted  as  a  function  of  dimensionless 
electron  energy,  E/EQ,  for  values  of  U  =  0,  1,  3,  5,  10, 
20.  Higher  values  of  U  correspond  to  an  increased  step¬ 
like  structure. 


4.  ONE-DIMENSIONAL  DENSITY  OF  STATES 

Since  the  Hamiltonian  in  equation  (1)  conserves  transverse  electron  momen¬ 
tum,  k  and  k  are  good  quantum  numbers,  and  an  electron  placed  in  a  state  of 
a  y 

definite  k^  will  remain  in  this  state  indef  ini  tely .  With  this  in  mind  we 
consider  the  Green's  function 

c,  E)  -  l  l  - ^~r  .  (30) 

a  k  E  -  E^  MO 
a  k 


in  which  the  summation  on  k^  is  omitted.  Using  this  Green's  function  we 
define  the  1-d  local  DOS  (including  both  spins)  as 


D^-d[z'Ezi  =  '  *  Im  Gi-d^r’r:E] 


1  /  ■h2kf  tf2k* 

—  I  l  Uak  (z)|2  6\E - — ^ 

L2  a  ka  a  a  \  2mQ  2mQ 


2mc  2mc 


where  ip  is  given  in  equations  (11)  and  (12),  and  E„  -  E  -fi2k|/2m  . 


The  function  D-|_d(z,E2]  gives  the  number  of  states  (labelled  by  a  and  kQ) 


per  unit  volume,  for  a  given  k^.  This  function  displays  peaks  which  are 


associated  with  the  resonances.  Substituting  the  explicit  form  of  the  wave- 


functions  into  equation  (31),  we  find 


Vd^2’^ 


2aL2E  /c 
o 


[Fe(§  ^)cos2(g  /e)  +  Fo(f  ^)sin2(|f  /t)]  , 


where  the  dimensionless  z-component  of  energy,  e,  is  defined  by  e  =  EZ/EQ. 
Again,  rather  than  looking  at  this  in  more  detail,  we  consider  the  integral  of 
D1_d(z,Ez)  over  the  well  volume 


■  1  D.-d<z’E2) d: 


2E  /e 
o 


-  f.d  ^ 


>in(  ir/e ) 


/it  /—\ 

>\2 


This  function  specifies  the  number  of  states  in  the  well  labelled  by  a,  k  , 


per  unit  energy,  for  a  given  k^.  In  the  limit  of  weak  barriers,  U  -*■  0 ,  the 


function  N-j  _d  C  Ez )  -»  1/(E0el/2)  which  is  the  DOS  for  a  1-d  free-electron  gas. 
In  the  limit  of  strong  barriers,  U  ■*  00 ,  the  number  of  states  in  the  well  per 
unit  energy  is  just  a  sum  of  delta  functions.  In  this  limit  the  resonant 
states  are  the  eigenstates  of  the  1-d  particle  in  a  box  problem,  and  the 
resonance  peaks  shift  to  the  appropriate  limiting  eigenvalues.  For  inter¬ 
mediate  values  of  U  (and  energies  0  <  c  <  U)  the  function  N1_d(Ez)  is  approxi¬ 
mately  a  sum  of  Lorentzians  (see  fig.  3).  The  lowest  energy  peak  is  composed 
predominantly  of  even  wavef unct ions ,  while  the  second  peak  is  composed  mostly 
of  the  odd  wavef  unctions .  In  figure  4  we  show  the  DOS  N-|_d(E)  for  several 
values  of  U,  in  the  energy  region  of  the  lowest  resonant  level.  The  inset  of 
figure  3  shows  the  peak  position  of  the  lowest  resonance  as  a  function  of 


>  V  " 


ft  rrr 


barrier  strength  U.  For  energies  much  larger  than  U,  c  >>  U,  the  f 
F  (ir/e/2)  -*•  1  and  the  DOS  returns  to  its  value  in  the  absence  of  b. 

M  ^  f  r  i  j.  i/fr  ./Cl 


Figure  4.  One-dimensional  density  of  states,  2  EQN^a(Ez) , 
plotted  versus  EZ/EQ  in  region  of  resonance  peak  for 
several  values  of  U. 


SUMMARY 


Within  the  context  of  a  simple  model  for  a  double-barrier  structure,  we 
solved  for  the  normalized  eigenstates.  Using  these  eigenstates  we  calculated 
the  3-d  local  DOS  between  the  barriers.  This  quantity  shows  a  crossover  from 
a  3-d  square  root  of  energy  behavior  to  a  quasi-2-d  staircase-like  behavior, 
as  the  barrier  strength  U  is  increased.  For  electron  energies  e  >>  U  the  DOS 
always  returns  to  the  free  electron  DOS.  We  also  calculated  the  1-d  DOS  for  a 

-f 

given  transverse  momentum  k^.  This  quantity  shows  sharp  peaks  at  energies 
corresponding  to  the  resonant  states.  In  a  more  realistic  model,  one  can  use 
the  width  of  the  lowest  peak  in  N-|_d(Ez)  to  reliably  estimate  the  lifetime  of 
the  lowest  quasi-bound  state.  The  inverse  of  this  lifetime  gives  an  estimate 
of  the  characteristic  frequency  above  which  the  resonant  contribution  to  the 
current  becomes  negligible. 
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